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FOREWORD 


This report covers the work completed on the research project 
"Chemoviscosity Modeling for Thermosetting Resins" for the period ended 
January 31, 1987. The work was supported by the NASA Langley Research 
Center (Polymeric Materials Branch of the Materials Division) under research 
grant NAG-1-569. The grant was monitored by Mr. Robert M. Baucom of the 
Polymeric Materials Branch. 
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STUDIES ON CHEMOVISCOSITY MODELING FOR THERMOSETTING RESINS 


By 

J. M. Bai 1 , T. H. Hou 2 , and S. N. Tiwari 3 

SUMMARY 

A new analytical model for simulating chemoviscosity of thermosetting 
resins has been formulated. The model is developed by modifying the well- 
established Williams-Landel -Ferry (WLF) theory in polymer rheology for 
thermoplastic materials. By introducing a relationship between the glass 
transition temperature Tg(t) and the degree of cure a(t) of the resin system 
under cure, the WLF theory can be modified to account for the factor of 
reaction time. Temperature-dependent functions of the modified WLF theory 
constants C^(t) and C£(T) were determined from the isothermal cure data. 
Theoretical predictions of the model for the resin under dynamic heating 
cure cycles were shown to compare favorably with the experimental data. 

This work represents progress toward establishing a chemoviscosity model 
which is capable of not only describing viscosity profiles accurately under 
various cure cycles, but also correlating viscosity data to the changes of 
physical properties associated with the structural transformation of the 
thermosetting resin systems during cure. 
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Mechanics, Old Dominion University, Norfolk, Virginia 23508. 

2 Adjunct Assistant Professor, Department of Mechanical Engineering and 
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Chapter 1 
INTRODUCTION 

Perhaps the single most Important property of a polymer with regard 
to specifying its processing characteristics is its viscosity. Viscos- 
ity governs the resin flow characteristics. For thermoplastic 
materials, the viscosity is influenced by local flow geometry and can 
vary with processing temperature and shear rate. Viscosity control 
becomes more critical and difficult in the processing of thermosetting 
resins, because of the onset of chemical reactions and the generation of 
heat during cure which causes the viscosity to vary with time. The term 
chemoviscosi ty refers to the variation of viscosity due to the polymer 
chemical reaction. The study of chemoviscosi ty is generally called 
chemorheology [1]*. 

In a typical autoclave operation for the fabrication of composite 
materials, the viscosity- time profile must allow first for the bubbling 
off of trapped gas, then fiber compaction with resin flow, and finally, 
laminate consolidation under applied pressure before the resin has 
gelled and ceases to flow. A typical cure cycle for the MACNAMITE 
AS4/3501-6 graphite pregreg (an amine-cured epoxy resin, Hercules 3501- 
6, reinforced with unidirectional graphite fibers [2]) is shown in Fig. 
1 . 1 . 


*The numbers in brackets indicate references. 


1 



500 


2 
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Fig. 1.1 A typical cure cycle for the MACNAMITE AS4/3501-6 graphite pregreg 
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At the start of a cure cycle, the polymerization reaction begins 
and the polymers chains lengthen. The increase in viscosity of the 

resin due to polymerization is, however, largely offset by the increase 
in temperature which is introduced through the cure cycle and the heat 
of reaction. Consequently, a drop in resin viscosity of two or three 
orders of magnitude at the initial stage of cure is not uncommon. As 
the resin continues to polymerize, crosslinks are formed and the 
viscosity of the resin system starts to increase at a faster rate. 
Finally, the rate of increase of viscosity approaches infinity at 
gelation. Therefore, any chemo viscosity profile is closely related to 
the reaction kinetics of the resin systems and the cure cycle 
(temperature profile) during processing. A fundamental understanding of 
the chemoviscosity kinetics of thermosetting resin has been necessary in 
order to obtain a precise processing control. 

The objective of present research is to establish a chemoviscosity 
model which is capable of not only describing viscosity rise profiles 
accurately under various cure cycles, but also correlating viscosity 
data to the changes of physical properties associated with the 
structural transformation of the thermosetting resin systems during 
cure. 

The thermal analysis and rheological analysis on net resin system 
were conducted by the differential scanning calorimetry and rheometry, 
respectively, during the course of the study. By introducing a well- 
established viscosity-temperature relationship already existing in 
polymer rheology based on the information of thermal analysis and 
rheological analysis, an analytical model for chemoviscosity profile of 
thermosetting resin during cure was formulated. The viscosity 
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calculated by the proposed model could be correlated with the physical 
properties through the two temperature-dependent functions C^(T) and 
C 2 ( T) as well as the glass transition temperature. In the following 
chapters, the experimental measurements on the thermal and rehological 
properties are first introduced and the modeling approaches and results 
are then presented. 


Chapter 2 


LITERATURE SURVEY 

Numerous researchers have devoted considerable efforts to establish 
an analytical model for the chemoviscosity growth profile of thermo- 
setting resins under curing. There are two different basic approaches 
in chemoviscosity modeling of thermosetting resin. The most common 
approach to the problem has been by empirically specifying the viscosity 
linearly as an exponential function of reaction [3,4], or with the aid 
of DSC thermal analysis, the degree of cure [5]. 

The empirical model proposed by Roller [3] to correlate isothermal 
viscosity data for a curing B-staged epoxy resin can be expressed as 

An ri » An n Q + kt (2.1) 

where t) is the time dependent viscosity, n Q is the zero-time viscosity, 
k is an apparent kinetic factor, and t is the time. Taking the 

temperature dependence into account, n Q could be expressed as 

n = tj exp( AE/RT) (2.2) 

o ® 

where n is the pre-exponential constant, A£ is Arrhenius activation 

oo T) 

energy of viscous flow, and R, the universal gas constant. Similarly, 
one can also express k with an Arrhenius expression 

k * k exp(-A£. /RT) (2.3) 

« k 


5 


6 


Here AE^ is the Arrhenius activation energy of chemical reaction. 

Substituting Eqs. (2.2) and (2.3) into Eq. (2.1), we have 

AE AE, 

Y% I/ 

An n(t) = An n + <— + t k exp(- -?-) (2.4) 

In the dynamic heating case, cure cycle is a function of time, i.e., T = 

T(t), and therefore, Eq. (2.4) can be generalized to 

AE t AE k 

An r)( t,T) = An t) + —2 + / k exp(- dt (2.5) 

00 K I 00 K I 

0 

In comparing the calculated viscosity by Eq. (2.5) to the 
experimental values. Roller [3] noted that there exists a large 
discrepancy, especially in the higher heating rate cases. It seemed 
that such empirical simulation of chemoviscosity is inadequate in 
describing accurately the non-linearity of the chemoviscosity-time 
profile associated with the advancement of resin under dynamic heating 
cure cycle. The model is also found to be very much batch-specific and 
cannot be related directly to the resin chemistry. 

In a similar manner. Carpenter [4] adopted the model given in Eq. 
(2.1) to simulate chemoviscosity of Hercules 3501-6 net resin system, 
which was cured under various conditions. The kinetics factor k in Eq. 
(2.1) was expressed as a different function of temperature than that of 
Roller. When comparing with the experimental values, however, similar 
inadequancies as those discussed above on Roller's modeling were also 
noted. 

Efforts were also made to model the time-dependence viscosity by 
Lee et al. [5]. The degree of cure of the reactive resin system during 
cure was introduced in chemoviscosity modeling by means of DSC. The 
model for isothermal case has been proposed as follows: 
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An t) * A + k a( t) (2.6) 

where AE 

A = In ^ 

AE^ is the activation energy for viscous flow and k is a constant which 
is independent of temperature. The viscosity in the dynamic case was 
obtained from Eq. (2.6). The temperature was calculated from the 
heating rate at every 0.5 seconds and the values were determined by 
numerical integration of reaction kinetics which was determined as a 
priority in the isothermal case. When compared with experimental data, 
it is seen that Eq. (2.6) is capable of describing the isothermal 
viscosity well; however, it is inadequate to model the non-linearity of 
chemoviscosity under the dynamic cure case. This approach, however, has 
suggested a way of relating chemoviscosity to a physical property, i.e. 
the degree of cure, which represents the various structure status of a 
reactive resin system under cure. 

The second approach is based on a modification to the well- 
established viscosity-temperature relationship existing in polymer 
rheology for thermoplastic materials. The parameters in such an 
equation can be expressed in terms of polymerization kinetics, and the 
chemoviscosity profiles as a function of reaction time can then be 
modeled for thermosetting resin system. 

The viscosity- temperature relationship, which is based on the 
Will iams-tandel -Ferry (WLF) theory [6], is given by 


n - C (T-Tg) 

Log ( J. ) = f 

% 


C 2 + T-Tg 


(2.7) 
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where Cj and C 2 are two material constants, and n-pg represents the 

viscosity at the glass transition temperature Tg of the given 

thermoplastic material. The normal use of the WLF equation for 

thermoplastic materials requires that Tg be constant while the 
temperature T is varied within 100 K above Tg for the specific polymer 
under study, and Tg can be replaced by any reference temperature Ts 

within the experimental range. It has been extensively documented in 
the literature that temperature-dependent viscosity of thermoplastic 
materials can be accurately described by the WLF equation within 100 K 
above its glass transition temperature Tg. However, during the curing 
of thermosetting resins, the monomers are initially polymerized and 
crosslinks are formed later. This is a system where Tg ( t) is changing 
and the curing temperature T is held constant (in an isothermal cure 
case, for example). The glass transition temperature Tg rises 
continuously and may eventually approach the curing temperature. Over 
the entire curing cycle, the material structure actually undergoes 
continuous phase transformation from a low molecular weight liquid to a 
high molecular weight polymeric melt, and eventually transforms to a 
crosslinked network. For a thermosetting resin system under cure, it is 
reasonable to assume that Tg ( t) of the material is always within 100 K 
lower than the cure temperature. The WLF equation should then be 
applicable to every Instant state of the resin system under curing. The 
central theme in modifying the WLF theory to describe the chemoviscosi ty 
of thermosetting resin is to take into account the reaction time factor. 

The applicability of the modified WLF theory in chemoviscosi ty 
modeling for thermosetting resin has been studied by Tajima and Crozier 
[1,7], Apicella et al. [8] and Hou [9] among others. 
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Tajima and Crozier [1,7] modeled the chemoviscosity as a function 
of polymerization kinetics for an amine (DDS)-cured epoxy resin by a 
modified WLF equation. The advancement of reaction was determined by 
the concentration of the hardener DOS. The parameters Ts and r)(Ts) in 
Eq. (2.7) were allowed to vary with times. The two material parameters 
Cj and C 2 were assumed to be independent of reaction time. A fairly 
good agreement between calculated and experimental data was observed. 

In his investigation, Hou [9] proposed an analytical model which 
was also based on modification of the WLF theory. Major assumptions 
involved were that the rate of reaction is diffusion controlled and was 
linearly inversely proportional to the viscosity of the medium over the 
entire cure cycle. The modified WLF equation became a first order 
ordinary nonlinear differential equation. Numerical solutions were also 
shown to compare favorably with the experimental results for several 
thermosetting systems under isothermal and dynamic heating cure 
conditions. It was demonstrated [9] that by such modified WLF equation, 
an analytical model with high degree of accuracy for the chemoviscosity 
simulation for any thermosetting resin system under various cure cycles 
can be established. The physical significances of the material 
parameters selected for the model were, however, difficult to extract 
for the particular resin system under investigation. 

A new approach in chemoviscosity modeling used in the present study 
will be formulated and discussed later in Chap. 5. 



Chapter 3 


EXPERIMENTAL METHODS 

3.1 Thermal Analysis 


3.1.1 Sample Preparation 

Sample ranging from 3 to 7 mg in weight was encapsulated in a 
standard aluminum sample pan. To prevent the sample holder being 
contaminated by the splitted sample, a small hole was made on the center 
of the lid. Tweezers were used throughout all the time for sample 

handling, loading and storing. The weight of the sample was measured 
before and after the test. The weight losses during the tests were 

found to be negligible. The prepared samples were stored in a 

refrigerator and removed just prior to the measurement. 

3.1.2 Appara tus 

Perkin-Elmer differential scanning calorimeter (Model DSC-2) 
coupled to a high quality po ten ti one trie recorder (Model 3314, a product 
of SOLTEC CORP.) was used in thermal analysis. 

The DSC is a sophisticated instrument for the measurement and 

characterization of the thermal properties of materials. It consists of 
an analyzer desk, draft shield, sample holder cover kit, concealed 

control panel cover, vacuum pickup device, etc. 
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When a transition such as melting, boiling and dehydration of 
crystallization occurs in the sample material, an endothermic or 
exothermic reaction takes place. The change in power required to 
maintain the sample holder at the same temperature as the reference 
temperature (i.e. its programmed temperature) during the transition is 
recorded as a peak. The chart abscissa indicates the transition 
temperature and the peak area indicates the total energy transfer to or 
from the sample. This direct calorimetric measuring principle of the 
instrument requires that each sample holder has a built-in heater and a 
temperature sensor. 

In the present work, the DSC-2 was upgraded to provide extra 
cooling capacity which can support operations down to -40°C by means of 
Intercooler II. In order to prevent moisture condensations on the 
sample holders while operated In the subambient environment, a package 
of intermediate range subambient accessory has also been installed. The 
heating chamber is enclosed in a dry box, and the sample holders can 
only be accessed through dry box gloves. A nitrogen purge gas line has 
been connected as well. Dry box is always purged with nitrogen before 
measurements started. The apparatus was Installed in an air-conditioned 
room maintained at 20 # C and 502 relative humidity. 

3.1.3 Cali bra tion 

The DSC-2 calibrations include the following: temperature 

(abscissa), energy (ordinate), a conversion constant K determination and 
the baseline optimization. All calibrations were performed by following 
the procedure recommended by the manufacturer [10]. 
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Three standard materials. Indium, Zinc and Potassium Chromate were 
used in the calibration. The characteristics of these three materials 
are given as follows: 


Transition Point (K) Transition Energy (cal/gm) 


Indium 

429.78 

6.80 

Zinc 

692.65 

25.90 

K2Cr04 

943.7*0.5 

8.50 


The temperature calibration is accomplished with the TEMP CALIB and 
TEMP CALIB RANGE controls of DSC-2. The construction for determining 
transition temperature is shown in Fig. 3.1. 


By scanning standard materials dynamically at a given rate, 
transition can be observed and the transition temperature can be 
determined. The comparison between measured temperature and actual 
reference temperature can then be made. The following equation was used 
to determine the change required in the setting of the TEMP CALIB RANGE 
control . 


where 

AR 


AR = Rr 


2 


r AT IND ’ AT ACT n r 1 

L AT J L 1000 J 

ACT 


(3.1) 


number of divisions change required to adjust TEMP 
CALIB RANGE control 


Rr setting of TEMP CALIB RANGE control 

ATj N 0 difference between measured transition temperatures 

AT A £j difference between the actual reference temperatures 
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By means of TEMP CALIB RANGE control adjustment, the difference 
between aTj N q and AT^y can be minimized. The calibration on 

temperature scale is completed when (ATj^p - ^ACT^ 1S ^ ess ^ an * 0*5 
K. The observed transition temperature for all the standards should be 
within 2 K of their actual values. A typical thermogram of K2Cr04 and 
Indium is shown in Fig. 3.2. In the present work, the deviation of 
measured temperature and actual temperature was within ± 1 K, which is 
quite satisfactory for our measurement. 


The energy calibration involves the determination of calibration 
constant K and is accomplished with the AW BALANCE control. In 
practice, the measurement of energy with the Model DSC-2 will 
necessarily involve considerations of an instrument calibration 
constant, the recorder chart speed, the sensitivity used, and the unit 
employed for area measurement, etc. To measure the total energy 
associated with a transition or reaction, it is necessary to integrate 
over the area under the peak with respect to time 




(3.2) 


The area of a peak is directly proportional to the energy of the change 
per unit weight, AH, in calories per gram, the sample weight, W and 
the chart speed S. It is also inversely proportional to the instrument 
RANGE setting, Rr. Therefore, we have 


„ 1 AHxWxS 

A IT — FT" — 


(3.3) 


where K is a proportionality constant and R' is setting RANGE of DSC. 
From the calibration with a standard material like Indium or Zinc 
through the adjustment of the AW calibration control, one can obtain 


K2Cr04 H h Indium 
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Fig. 3.2 The DSC thermograms of K2Cr04 and Indium under dynamic heating 
condi tion. 
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the value of K from Eq. (3.3). Table 3.1 shows the results obtained for 
calibration constant K used in the present work. 

The baseline optimization is accomplished with SLOPE control and 
AT BALANCE control. By adjusting AT BALANCE and SLOPE control, the 
curvature and the decline of baseline could be improved until the 
satisfactory flatness of the baseline is obtained. Since the baseline 
affects the area of thermogram measured, it is recommended that the 
baseline optimization be conducted for different scanning conditions and 
materials. 


3.1.4 Measurements of Degree of Cure a(t) 

The degree of cure was defined as the fraction of heat released, up 
to time t for the resin system under cure [5] measured by Perkin-Elmer 
DSC-2. The measurement was made under the dynamic cure condition with a 
given heating rate. During the temperature scan, the rate of heat 
generated was recorded and plotted versus time and the area under the 
curve provided the amount of heat released (Fig. 3.3). The total heat 
of reaction, Hy, is calculated by the expression 




(3.4) 


where ty is the time required to complete the reaction, 
released up to time t, H, could also be calculated as 

U dH Ht 

H = J . dt 

z o 

and the degree of cure can be determined by 


The heat 


(3.5) 
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Table 3.1 DSC calibration constant determined from 
the thermogram of standard material 
Indium (weight = 5.0 mg, AH = 6.8 cal/gm 
and heating rate = 20 K/min) v 


Chart. Rec. 
Speed 
sec/in. 

Range 
meal /sec 

Area 

in^ 

K 

1/in 

5 

10 

7.0 

0.97 

5 

20 

3.53 

0.96 

10 

10 

3.55 

0.96 

20 

10 

1.76 

0.97 



i 


19 


a 



(3.6) 


The rate or degree of cure can also be calculated by the expression 


do _ 1_ dH 
dt ‘ H t dt 


(3.7) 


The measurements of degree of cure were made on Hercules 3501-6 net 
resin system under dynamic cure conditions according to the procedure 
outlined above. The heating rates were 10, 20, 40 and 80 K/min 
respectively. The total heat of reaction, H, is calculated by a 
numerical integration scheme based on Simpson's rule [111. Experimental 
results are presented in Chap. 4. 


3.1.5 Measurements of Glass Transition Temperature Tg(t) 

Measurements of glass transition temperature under both dynamic and 
isothermal condition were performed. Dynamic curing experiment was 

performed as follows: The sample was first heated at heating rate of 20 

K/min from 260 K to a specified temperature. A rapid quench to the 
initial temperature then followed. The resultant thermogram was 

obtained by re-scanning the sample at a rate of 20 K/min. A new desired 
temperature could be reached in this stage and subsequent measurements 
could be performed on the same sample. The experimental cure cycle is 
shown in Fig. 3.4. 

In the present study, glass transition temperature at eight 
different states (during the advancement of the reactive resin system 
cured at a constant rate of heating condition of 20 K/min) were 
measured. The measurements were repeated for ten different samples. 


9 } 



Fig. 3.4 The experimental cure cycle of glass transition temperature measurements 
for Hercules 3501-6 under dynamic heating condition. 
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In the Isothermal scanning measurement, the experiment was 
performed as follows: The sample was first heated at a rate of 320 

K/min from 260 K to a specified curing temperature. At the completion 
of various holding times at the temperature, a rapid quench to initial 
temperature followed. Subsequent thermograms were obtained by re- 
scanning the quenched samples at a rate of 20 K/min from 260 K to a 
temperature above its Tg. A number of measurements of glass transition 
temperature under isothermal scanning condition was conducted. The 
selected curing temperatures and the length of curing times were 
determined according to the reaction kinetics reported by Lee et al. 
[5]. 


3.2 Rheological Measurements 


3.2.1 Sample Preparation 

All samples were preconditioned at room temperature by pressing 
about 30 gm of resin between two pieces of non-porous Teflon cloth (CHR- 
6TB) to produce a resin wafer with a thickness of 1.25 mm. The wafers 
were stored in a vacuum oven for a couple of hours at a temperature of 
about 340 K until the volatile bubbles were no longer evident. The 
sample preconditioning was used to remove moisture and any other 
volatiles that might create erratic variations in rheology. After 
degasing, the sample was stored in the refrigerator and removed just 
prior to the experimental measurement. A disk sample was prepared by a 
circular puncher with a diameter of 2.50 cm. 
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3.2.2 Appara tus 

The Rheometrics System Four [12] was used to measure the 
rheological properties of the material. 

The Rheometrics System Four is a fully automated laboratory 

instrument for characterizing the rheological properties of a broad 
range of materials including solids, fluids and melts. It is a modular 
designed system of four independent servos linked to a common computer 
system for control and data acquisitions. These four servo motors are: 
Dynamic, Linear, Fluids and Steady which are selected for optimum 

response in their respective modes allowing a broad range of test 
geometries and deformation histories. 

The dynamic servo motor and a parallel plate test geometry are 

selected in the present study to measure chemo viscosity of the reaction 
epoxy resin system. Sample temperature is controlled with a convected 
gas environmental chamber. The torque and normal force generated in 
response to the imposed motion is measured by the transducer. The 

computer calculates stresses from the torque and normal force 
measurements and combines these values with measured sample deformation 
to calculate values of the selected rheological properties. The 

computer and associated electronics also provide precision control and 
rapid data analysis with results conveniently printed out in graphical 
and tabular forms. 

In measuring chemoviscosi ty of Hercules 3501-6 epoxy resin, a time 
sweep mode or cure mode were selected respectively in isothermal or 
dynamic cure condition. 
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3.2.3 Calf bra tion 

The calibration of Rheometrics includes: transducer normal and 

torque calibration, temperature calibration and phase calibration. 

Transducer calibration must be performed whenever a transducer is 
changed. Periodic calibration checks should be performed to ensure the 
accuracy of test results. Calibration is accomplished through adjust- 
ment of normal calibration control and torque control by comparing 
measured data with standard weight. 

The phase calibration is accomplished to compensate for phase shift 
in the electronics through the adjustment of appropriate potentiometer. 
Calibration was performed with the dynamic motor in TORSION RECT using 
an elastic material such as steel which possesses a phase angle of zero 
degrees and a low viscosity newtonian fluid which possesses a phase 
angle of -90 degrees. Calibration in TORSION RECT or BENDING with a 

material such as steel is the preferred method as the phase angle of 
steel is not frequency or strain dependent. For newtonian fluids and 
elastic samples, the tolerance is -0.4 and -0.25 degrees respectively. 

3.2.4 Measurement of Chemoviscosity 

The chemoviscosity data under isothermal condition has been 
measured previously [13]. The selected temperatures were 360, 375, 385, 
399, 410, 425 and 435 K. 

The measurement of chemoviscosity in a dynamic case was conducted 
in present study. During the measurement, the sample was confined to 
the gap between two parallel plates mounted in the Rheometer. The top 
plate was motor driven about its axis while the bottom plate was mounted 
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on a torque transducer for force measurement. The typical gap between 
the parallel plates was 1.2 mm. The plates and sample were enclosed in 
a heating chamber where temperature control was provided. Selection of 
the strain value was to assure that the measurements were performed 
within the material's linear viscoelastic response range. 

It is well known that the viscoelastic properties of polymeric 
materials respond differently to different frequencies, and selection of 
the frequency value was made to assure that the complex viscosities 
measured were within the Newtonian region [12]. The cure mode of the 
Rheometer was used during testing which automatically programmed 
different temperature profiles as the epoxy was cured. 

The sample was first kept at 300 K for about 5 minutes before 
measurement started. Experimental measurements were performed under 
different rates of heating conditions. 



Chapter 4 


EXPERIMENTAL RESULTS 


4.1 Material 

The material used in the study is Hercules 3501-6 net epoxy 

resin. This thermosetting system was selected because of its wide use 
in aerospace industry. The material received from the manufacturer was 
used without further treatment, and was kept in the dark and cold 
conditions at around -5*C. 

The material contains the following standard concentrations of 

starting components [14]: 


COMPONENT TOTAL PERCENT (wt) 


Tegraglycidyl Methylenedianiline 

56.5 

Alicyclic Diepoxy Carboxylate 

9.0 

Epoxy Cresol Novalac 

8.5 

4,4* Diaminodiphenyl Sulfone 

25.0 

Boron Trifluoride Amine Complex 

1.1 


4.2 Experimental Data of Thermal Analysis 
4.2.1 Degree of Cure a(t) 

4. 2. 1.1 Isothermal Case . Lee et al. [5] performed thermal 

analyses on the same epoxy resin system by DSC. The degree of cure a 
and the rate of change of degree of cure da/dt were determined from 
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i 

I 

I 

the results of the isothermal scanning experiments. Eight different 
temperatures were selected, and the results are reproduced in Fig. 

4.1. It is noted from the figure that both higher values of the degree 

I 

| of cure and rate of change of degree of cure are decreasing functions of 

cure time t, and eventually level off at a low value (< 0.1) when a high 
degree of cure is reached. The following equations were found to 
describe the temperature-dependent data rather accurately: 

da/dt = (K +K a) (1-a) (B-a) a < 0.3 (4.1a) 

n l n 2 

i 

I 

j da/dt * K (1-a) > 0.3 (4.1b) 

n 3 

where 

I K n = A i exp(-AEl/RT) 

K * A, exp(-AE2/RT) 

i n 2 <L 

K = A. exp(-AE3/RT) 
n 3 J 

The values of constants in Eqs. (4.1a) and (4.1b) are summarized in 
Table 4.1. Comparisons between the measured values of da/dt and the 
values calculated by Eqs. (4.1a) and (4.1b) are also shown in Fig. 4.1 
by symbols and solid lines respectively. 

4. 2. 1.2 Dynamic Case . The measurements of heat of reaction and 
degree of cure under dynamic heating condition were performed with four 
different heating rates, 10, 20, 40 and 80 K/min respectively. A 
typical thermogram at 20 K/min is shown in Fig. 4.2. It is seen that 
there are two "humps" in the curve. These humps are caused by two 
different major reactions occurring during the cure. The total heat of 


RATE OF DEGREE OF CURE, da/d I 


Colcutoted 


425 K 


430K 



440K i£> 


[•]! 


465K 


470 K 


0.4 0.6 0 0.2 0.4 06 

DEGREE OF CURE, a 


Fig. 4.1 Rate of degree of curve vs. degree of cure reported by Lee 
Loos and Springer [5] for Hercules 3501-6 resin system 
under isothermal cure cycles; solid lines were calculated 
using Eqs. (4.1a) and (4.1b). 



Table 4.1 Yalues of constants Eqs. (4.1a) and (4.1b) for 
Hercules 3501-6 resin system under isothermal 
condition 


B = 0.47 

Aj * 2.101 x 10 9 min"* 
A 2 = -2.014 x 10 9 min" 1 
A 3 = 1.960 x 10 5 min" 1 
AE^ = 8.07 x 10 4 J/mole 
AE 2 = 7.78 x 10 4 J/mole 
AE 3 = 5.66 x 10 4 J/mole 


Temperature 



opug 


oxg 


Fig. 4.2 A typical rate of heat generation vs. time curve of Hercules 3501-6 
resin system provided by dynamic scanning experiments under 20 K/min 
rate of heating condition. 
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reaction H was calculated by numerical integration according to the 
Simpson's rule. The algorithm was programmed in FORTRAN and was run by 
a VAX computer system. Program codings are included in Appendix B for 
references. Calculated results are tabulated in Table 4.2. It can be 
seen that Hy is independent of heating rates. An average value of H t = 
120±5.0 cal/gm was chosen to be used in all future calculations. The 
total heat of reaction calculated by Lee et al. [5] in dynamic scanning 
experiment was 113.4 cal/gm. A ten percent difference between the two 
results is noted. This is possibly caused by materials with different 
aging histories and the scheme employed in determining thermogram and 
the area under the peak of the curve. 

The degree of cure o(t) was also calculated based on the 

thermograms. The results are plotted in Fig. 4.3. Open and filled 
symbols represent results from different runs under the same condition. 
The reproducibili ty of the measurements appears to be very satisfactory. 
Numerical solution of Eqs. (4.1a) and (4.1b) under dynamic heating 
condition where temperature was not kept constant is outlined as 
follows. 

In a dynamic heating case, reacting kinetics of resin system 
becomes a first-order, non-linear and initial-value differential 
equation, i.e., 

|f • fit, a] . a(t 0 ) • a 0 (4.2) 

Assume that f is sufficiently differentiable with respect to either t or 
a. It is known [15] that Eq. (4.2) possesses a unique solution if 
is continuous on the interval of interest. If a(t) is the exact 

solution of Eq. (4.2), then we can expand a(t) in a Taylor series 

about the point t * t 0 as 
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Table 4.2 Total heat of reaction for Hercules 3501-6 
resin system measured at different rate of 
heating conditions 


Weight 

mg 

Rate 

°K/min. 

Hea ting 
Speed 
sec/in. 

Chart. Rec. 
Range 
meal /sec 

h t 

cal/gm 

4.6 

10 

60 

0.5 

118.702 

5.5 

10 

60 

0.5 

120.129 

6.4 

10 

60 

0.5 

126.472 

5.2 

20 

40 

0.5 

119.474 

6.7 

20 

40 

0.5 

119.474 

5.7 

40 

20 

1.0 

125.37 

7.6 

80 

10 

2.0 

118.34 

5.3 

80 

10 

2.0 

118.32 


Heating Rate (°K/min) 

Ave. H t 

(cal/gm) 

10 

121.77 

± 4.136 

20 

122.39 

± 4.131 

40 

120.02 

± 6.26 


80 


118.3 ± 2.21 



Fig. 4.3 Degree of cure a(t) measured for Hercules 3501-6 resin system cured under 
four different rates of heating condition shown. 
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a(t) 


“o + (t • 'o’ 


b «<V 

at 


(t - t Q ) 2 a 2 a(t Q ) 


7 T 


at 


7 


(4.3) 


By taking the total derivative of Eq. (4.2) with respect to t, (keeping 
in mind that f is an implicit function of a), we have 


da 

at 


= f(t, a) 


(4.4) 


a 2 a _ fl _ _3f , 3f 3a 
77 'at aa at 


(4.5) 


For a < 0.3: 


f(t, a) * (K + K a) (1 - a) (B - a) 

n l n 2 


ftf AE 1 &9 

9 (K + -nr- K a) (1 - a) (B - a) 


* ' is? ,R »r *i \ 


af 


^ = K (1 - a) (B - a) - (K + K a) (B - a) - (K + K a) 
oo n 3 n^ ^2 n l n 2 


For a > 0.3: 


_ AE 3 

f(t, a) = A 3 e ^ (1 


a) ■ K n (1 - a) 
n 3 


af . aE 3 n . 

-rr T K (1 - a) 

5t Rkt 2 n 3 


af 

3a 


= - K. 


By denoting a(t..) = a., we have 


2 2 

3a. At5 3 a. 

a i+l = a i + At i ~ + TT 77“ 

oz 


(4.6) 

(4.7) 

(1 - a) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 


(4.12) 
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a(t Q ) * a( 0) = 0 


where i - 0, 1, 2, 3,..., and 

5a . 


* 0^ + K a . ) (1 - a.j ) (B - ) a < 0.3 (4.13a) 


da. 

W ‘ K n„ {1 " a i } 


a > 0.3 (4.13b) 


*2 

6 a. 
5t 2 


AE, AE ? 

Wt (K n : + K n 2 V {1 


a j ) (B 


) + [K n (1 - a,. ) ( B - a.. ) 


- \ + % n> (B - V - \ * k „ 2 V (1 - V 3 1« 


+ K. o ( ) (1 - a.) (B - a ( )] 

i*2 ill 


a < 0.3 (4.14a) 


5 a- AE^ 9 

T~ = “7 K n (1 - a <) - K n (1 * “,*) 

dr Rkr n 3 1 n 3 1 


a > 0.3 (4.14b) 


The numerical results are plotted by solid curves in Fig. 4.3. 
Considering the fact that the temperature ranges of 400 to 475 K covered 
by isothermal experiments, as represented by Eqs. (4.1a) and (4.1b), is 
narrower than the ranges of 320 to 600 K covered here by the dynamic 
heating experiments, the agreements shown in the figure are rather 
satisfactorily. Figure 4.3 also implies that information such as 
chemoviscosi ty r|(t) and degree of cure a(t) under dynamic heating 

cure conditions can be related within certain accuracy to those obtained 
under isothermal curing conditions. 
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4.2.2 Glass Transition Temperature Tg(t) 

4. 2. 2.1 Isothermal Case . It was found that the major reaction of 
material occurred within the range of 420 to 500 K in dynamic curing 
condition. The range of degree of cure was between 0.1 and 0.5. It is 
the range in which the simulations of the chemoviscosity-time profiles 
are performed. 

The glass transition temperature Tg(t) was determined for 
isothermal curing conditions at five temperatures of 430, 440, 450, 460 
and 480 K respectively. A set of three degrees of cure was selected for 
each curing temperature as shown in Table 4.3. The corresponding set of 
curing time, t, was determined by solving Eqs. (4.1a) and (4.1b) 
isothermally, and this is presented by 


K 

a n 2 

t * An(l + |^“ ■ cc) 

n 2 n l 


- b An(l-a) - c An(l- -gO a < 0.3 (4.15a) 


where 


t = JL An(^) + tc 
n 3 


a = (B-l)/d 


b = (K + K B)/d 
n l n 2 


- (K + K )/d 
n l n 2 


a > 0.3 (4.15b) 


d = (K_ + K ) {K + (K 


n~ n. 


K ) B 
n 2 


- K } 


and tc is the time given by Eq. 4.4a for a - 0.3. The selected curing 
temperatures and the length of curing times for isothermal measurements 
are tabulated in Table 4.3. 
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A set of typical traces obtained from isothermal measurements is 
shown in Fig. 4.4. During the glass transition, the material undergoes 
a transition from glassy state to rubbery state, and the specific heat 
Cp of the material is changed. Thus, a step change is observed in the 
thermogram. Tg is determined at the arrow marks shown in the curve. 
The data obtained from the measurements is listed in Table 4.3. 

It is found that the glass transition temperature in isothermal 
curing condition depends not only on the degree of cure (or curing 
time), but also the curing temperature. A plot of Tg vs. a at various 
curing temperatures is given in Fig. 4.5. The data was fitted by linear 
Least Square technique. The relationship between glass transition 
temperature and degree of cure can be expressed as 

Tg = dl(T) a(t) + d2 (4.16) 

where dl is a temperature dependent parameter. 

An Arrhenius type plot of dl vs. 1/T is shown in Fig. 4.6. By 
means of linear Least Square fit, two straight lines are represented, 
respectively, by 

Logdl(T) = -639.56/T + 3.357 T < 450 K (4.17a) 

Logdl(T) = -3216. 49/T + 8.879 T > 450 K (4.17b) 

The correlating factors of the above two linear fits are 0.999. 

4. 2. 2. 2 Dynamic Case . The glass transition temperatures for the 
material at twelve different states during cure at a constant rate of 
heating condition of 20 K/min were measured. The measurements were 
repeated for ten different samples. The method of determination of Tg 


Hercules 3501-6 net resin 20 K/min. scan 
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opug 


IP/HP 


0X3 


Fig. 4.4 The typical thermogram of glass transition of Hercules 3501-6 resin system 
under isothermal cure condition; the curing temperatures are 430 and 450 K 
respectively. 



400 



Fig. 4.5 The glass transition temperature Tg as a function of degree of cure for 
Hercules 3501-6 resin system under different isothermal curing 
temperatures shown. 




m 



Fig. 4.6 Plots of temperature dependent parameter of Eq. (4.5) dl(T) vs. 1/T for 
Hercules 3501-6 resin system under isothermal curing conditions. 
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from the trace was the same as that in the isothermal case. The 
measured results are tabulated in Table 4.4 together with degree of cure 
a(t) measured by the procedures outlined above. Average values of 
Tg(t) and standard deviation at different resin states defined by curing 
time t (or curing temperature T) are included as well. It is noted that 
standard deviations increase for increasing scanning temperatures. 
Difficulties in determining Tg accurately come from two sources: (i) 

the higher the scanning temperatures, the less sharpness for the step 
change of dH/dt at the glass transition temperature and (ii) the 
existence of baseline slope which is usually distorted at higher 
temperature regimes. Nevertheless, a 90 % confidence level can be 
achieved with ± 3 % of the experimentally determined values of Tg ( t ) . 

Values of Tg ( t) are plotted in Fig. 4.7. For the dynamic case, 

cure cycle T(t) = 20t with a constant heating rate of 20 K/min is 
substituted into Eq. (4.3), the glass transition temperatures in dynamic 
heating condition can then be calculated from the formulas for the 
isothermal case. Comparisons between the calculated values (the solid 
curve) and experimental data are shown in Fig. 4.7. Favorable 
comparisons are seen for the changes of the glass transition temperature 
in dynamic heating case below 0.5 degree of cure, where the simulations 
of the chemoviscosi ty-time profiles are performed. 

4.3 Experimental Oata of Rheological Measurement 
4.3.1 Isothermal Case 

The chemo viscosity results obtained for isothermal conditions are 
shown in Fig. 4.8. The curing temperatures are 360, 375, 385, 399, 410, 
425 and 435 K. It is noted that for a higher cure temperature, a 



Table 4.4 Glass transition temperature measured during reaction advancement of 
Hercules 3501-6 resin system under 20 K/min rate of heating condition 
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Fig. 4.7 Plot of Tg(t) vs. a(t) for Hercules 3501-6 resin system under dynamic 
heating condition at 20 K/min. 



44 



(S0s.iocl) Xjisoosja 


Fig. 4.8 Plot of chemoviscosity vs. time t for Hercules 3501-6 resin system under 
seven isothermal curing temperatures shown. 
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shorter curing time is needed for reaching the same viscosity level. 
For a given time, t, the initial decrease of viscosity during the curing 
was attributed to the temperature effects. Apparently, the non- 
linearity of the isothermal viscosity vs. time exists for the material. 

4.3.2 Dynamic Case 

Experimental resul ts of chemoviscosi ty in the dynamic case are 
shown in Figs. 4.9-4.11. Three different dynamic heating rates of 2, 3 
and 5 K/min were used. The initial decrease in viscosity is caused by 
the increase in curing temperature. At the gelation, the viscosity is 

noted to increase exponentially. 


1E6 



Fig. 4.9 Plot of chemoviscosi ty n vs. time t for Hercules 3501-6 resin system under 
dynamic curing condition at 2 K/min. 
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Fig. 4.10 Plot of chemovi scosi ty n vs. time t for Hercules 3501-6 resin system 
under dynamic curing condition at 3 K/min. 



Hercules 3501-6 net resin 



rsesioch Aiisoosia 


Fig. 4.11 Plot of chemoviscosity n vs. time t for Hercules 3501-6 resin system 
under dynamic curing condition at 5 K/mln. 



Chapter 5 


CHEMOVISCOSITY MODELING 
5.1 Theory 

The Will iams-Landel -Ferry (WLF) theory [6] states that for Tg < T < 
(Tg+100 K), the temperature-dependent viscosity is given by the 

expression 


log (— ) 

% 


C^T-Tg) 
C 2 + T-Tg 


(5.1) 


where and C 2 are material constants. The temperature dependency of 
the viscosity as shown in Eq. (5.1) has been tested and found valid for 
materials ranging from dilute polymeric systems and thermoplastic melts 
to rubbers and elastomer. Equation (5.1) can be derived from the semi- 
empirical Doolittle equation which relates the viscosity to the free 
Volume of the liquid. This will be discussed later in this chapter. 

The normal use of the WLF equation for thermoplastic materials 
requires that the glass transition temperature Tg be constant while the 
temperature T is varied for the specific polymer under study. Equation 
(5.1) is applicable for a temperature T up to 100 K higher than the Tg 
of the material. However, during cure of thermosetting resin, the 
monomers are initially polymerized and crosslinks are formed later. 
This is a system where Tg(t) is changing and the curing temperature T is 
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held constant (e.g. the isothermal case). The glass transition 
temperature Tg rises continuously and eventually approaches the curing 
temperature. Over the entire curing cycle, the material structure 
actually undergoes a continuous phase transformation from the low 
molecular weight liquid to the high molecular weight polymeric melt, and 
eventually transforms to form crosslinking networks. If it is assumed 
that Tg(t) of the material is always lower than the cure temperature T, 
and that (T-Tg(t)) is always within 100 K, then the WIF theory should be 
applicable to all polymer structure phases during the cure. However, 
before Eq. (5.1) can be applied to describe the chemoviscosity of the 
curing resin, modifications have to be made to properly account for the 
reaction time, t, of the thermosetting resin system. 

One method of introducing reaction time factor into the WLF 
equation is to express the glass transition temperature Tg(t) as a 
function of a(t), the degree of cure. Several different techniques 
have been used in literature which include thermal, spec trop ho tome trie 
(FTIR) and chromatographic (HPLC) measurements to determine the extent 
of cure. In the present study, we assume that for the resin system 
under investigation, a(t) at time t is equal to the fraction of heat 
released, as measured by DSC, up to time t for the resin system under 
cure. The a's thus determined are a function of the curing 
temperature and time. The two material constants, Cj and C 2 , should be 
varied with reaction time because of the transformation of different 
material structures during cure. 
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5.2 Chemoviscosi ty Modeling 

The procedure used in modeling of chemoviscosi ties of Hercules 
3501-6 epoxy resin is described by the following two steps: 


1. Temperature-dependent material constants C^fT) and C 2 (T) in Eq. 
(5.1) are first determined by the experimental data for the 
material cured under isothermal condition. 

2. The chemoviscosi ty T)(t,T) established by isothermal cure data 
under procedure 1 is allowed to follow various dynamic heating 
cure cycles T(t). The viscosity values thus calculated are 
compared with the experimental data. 

5.2.1 Isothermal Case 

For the material cured isothermally, the relationship between 
degree of cure a(t) and curing time t is established by integrating 
the kinetics model of Lee et al. [5], As discussed in Chap. 4, this is 
given by 


a n 2 

t = |ji — An(l + ■ o) 

n 2 n l 


b An (1-a) - c An (1 - |) a < 0.3 (4.4a) 


t = “ TT“ * n ( ^T ) + tc 
n 3 


a > 0.3 


(4.4b) 


The glass transition temperature Tg(t) and the degree of cure a(t) is 
related experimentally by 


Tg = dl(T) a(t) + d2 


(4.5) 
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where 

Logdl(T) * -639.56/T + 3.357 T < 450 K (4.6a) 

Logdl(T) * -3216. 49/T + 8.879 T > 450 K (4.6b) 

d2 = 283.0 

Substituting Eqs. (4.4a), (4.4b) and (4.5) into Eq. (5.1), the 

calculated o(t) at a given isothermal cure condition can be made to 
fit the experimental data through the adjustments of the two material 
constants and C 2 . During the model calculation, Is assigned a 

value of 10*3 poises. This is the viscosity determined at Tg for many 
glass-forming substances and polymer systems [16]. 

The Cj(T) and C 2 fT) values determined for seven different 
temperatures are listed in Table 5.1, and the model predictions are 
drawn as a solid line in Figs. 5.1 to Fig. 5.7. It is noted that the 
model is capable of describing the non-linear relationship between 
viscosity which is below the order of 10^. This is considered a 

significant improvement over the empirical linear model used by Lee et 
al. [5]. 

An Arrhenius plot of Cj and C 2 respectively, vs. temperature, is 
shown in Figs. (5.8) and (5.9). The linear Least Square technique was 
used again and was found to fit the data reasonably well. The two 
straight lines are represented, respectively, by 

Cl = 4.067xl0 4 j - 54.252 


for T < 382 K (5.2a) 
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Table 5.1 Values of Cj and C£ determine for 
seven different temperatures under 
isothermal condition. 


Tc (K) 

C 1 

c 2 

360 

59.0 

300 

375 

53.5 

300 

385 

51.8 

300 

399 

85.7 

600 

410 

110.7 

850 

425 

128.7 

1100 

435 

133.4 

1300 



400 



Fig. 5.1 Chemoviscosi ty of Hercules 3501-6 resin system cured under isothermal 

curing temperature 360 K; triangles denote the experimental data; solid 
curve denotes the model prediction of Eq. (5.1). 
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Fig. 5.2 Chemoviscosity of Hercules 3501-6 resin system cured under isothermal 
curing temperature 375 K; triangles denote the experimental data; sol 
curve denotes the model predictions. 
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Fig. 5.3 Chemoviscosity of Hercules 3501-6 resin system cured under isothermal 

curing temperature 385 K; triangles denote the experimental data; solid 
curve denotes the model prediction of Eq. (5.1.). 
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Fig. 5.4 Chemoviscosity of Hercules 3501-6 resin system cured under isothermal 
curing temperature 399 K; triangles denote the experimental data; solid 
curve denotes the model prediction of Eq. (5.1). 



58 



Fig. 5.5 Chernov i sco si ty of Hercules 3501-6 resin system cured under isothermal 

curing temperature 410 K; triangles denote the experimental data; solid 
curve denotes the model prediction of Eq. (5.1). 
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Fig. 5.6 Chemoviscosity of Hercules 3501-6 resin system cured under isothermal 

curing temperature 425 K; triangles denote the experimental data; solid 
curve denotes the model prediction of Eq. (5.1). 
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Fig. 5.7 Chemoviscosity of Hercules 3501-6 resin system cured under isothermal 
curing temperature 435 K; triangles denote the experimental data; 
solid curve denotes the model prediction of Eq. (5.1). 
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C2 = 300 


Cl = 1 2.776xl0 5 j 779.09 


for T > 382 K (5.2b) 


C2 = -3.33xl0 6 y + 8959.65 

These temperature dependent Cj{T) and C 2 (T) are used to simulate the 
chemoviscosity-time profiles under dynamic heating cases. 


5.2.2 Dynamic Heating Case 

Similar techniques as those used in the isothermal cases are 
employed here. The advancement in the degree of cure under dynamic 
heating conditions can be calculated from Eqs. (4.1a) and (4.1b) by 
means of the Taylor's series expansion as discussed in Sec. 4.2. By 
substituting Eqs. (4.5), (5.2a) and (5.2b) into Eq. (5.1), the 

chemoviscosity-time profiles for the resin system cured under dynamic 
heating conditions can be simulated. 

Values of a(t) determined by the Taylor's series expansion depend 
slightly on the time interval At chosen. Theoretically, a(t) can be 
determined more accurately when a smaller At is used. The effect of 
step size At on the calculated degree of cure a(t) is given in Table 
5.2. The differences in a(t) are within 2% and are negligible. The 
step size of At*0.01 was chosen for all calculations. 

Chemoviscosities of the Hercules 3501-6 epoxy resin system under 
dynamic heating cure condition were obtained at heating rates 2, 3 and 5 
K/min, respectively. All cure cycle were started at 300 K. The results 
of the model predictions are shown as solid curves in Figs. 5.10-5.11. 


Table 5.2 The calculated degree of cure a(t) by 
different time step size At 


Curing Time (min) 


a(t) 


At=0.01 

At=0.05 

At=0.1 

6 

.150226 

.153258 

.157127 

8 

.347579 

.351781 

.357062 

10 

.606030 

.611808 

.619061 


12 


899822 


904560 


910432 


Hercules 3501-6 


net resin 


2 K/min. scan 


OOOOQ 


Time (min.) 


Fig. 5.10 Chemovi scosi ty of Hercules 3501-6 resin system under dynamic 
curing condition at heating rate of 2 K/min; circles denote 
the experimental data; solid curve denotes the model 
prediction of Eq. (5.1). 





Viscosity 


Hercules 3501-6 net resin 


3 K/min 


5 K/min 


1E2 


scan 


Time (min.) 


Time (min.) 


Chemoviscosity of Hercules 3501-6 resin system under dynamic 
heating condition at heating rate of 3 and 5 K/min. 
respectively; circles denote the experimental data; solid 
curves denote the model prediction of Eq. (5.1). 
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It is noted that the model calculations compared favorably with the 
experimental data for all three cases. 

The WLF model used in the present study was originally derived from 
the semiempirical Doolittle equation for the viscosity of a liquid as 
[17] 

V - V f 

Ann = AnA" + B" <— _I) (5.3) 

v f 

This gives an expression for the viscosity of a system in terms of two 
constants A and B. The quantity V is the total volume of the system 
while Vf is the free volume available to the system. Rearrangement of 
Eq. (5.3) gives 

Ann = AnA" + B" (| - 1) (5.4) 

where f is the fractional free volume Vf/V. It is assumed that above 
tne giass transition temperature, the traction tree volume increases 
linearly, that is 


f « fg + a f (T - Tg) 


T > Tg (5.5) 


where f is the fractional free volume at T, fg is the fractional free 
volume at Tg, and is the thermal coefficient of expansion of the 

fractional free volume above Tg. Substituting Eq. (5.5) into (5.3), the 
Doolittle equation becomes 



a 




(5.6) 


Equations (5.1) and (5.6) have the same form, and we can have 
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B' 1 

C1 = 2.303 Tg 


f 



(5.7) 


where B", a parameter in the Dolittle equation, is experimentally found 
to be close to unity. 


The concept of free volume is difficult to derive in a precise 
manner. In an approximate way, we can represent the segment of a 
polymer chain by rigid bodies and the free volume as the holes presented 
between these segments as a result of packing requirements. Presumably, 
the free volume reaches a constant value at Tg, which is too small to 
allow the large-scale conformational rearrangements of the chain 
backbones associated with Tg to occur. Above Tg, on the other hand, 
free volume increases and becomes sufficiently large to allow such 
motions to occur. 


The constants Cj and C 2 have been shown to vary from polymer to 
polymer [18]. During the curing of thermosetting resin, the mononers 
are initially polymerized and crosslinks are formed later. The material 
undergoes continuous phase changes from dilute polymeric solutions to 
polymer melts, and eventually the elastomers. A single set of Cj and C 2 
for each isothermal curing condition is found here, however, to describe 
the chemoviscosi ty-time profile adequately. This single set of values 
is therefore considered to describe a material state which represents an 
average over various structures occurring inside the material during 
cure. Values of and C 2 do vary with different isothermal curing 
temperatures as can be seen in Table 5.1. The fractional free volume 
and thermal expansion coefficient calculated by Eq. (5.7) are given in 
Table 5.3. 
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Table 5.3 The values of the fractional free volume and 
thermal expansion coefficient determined 
under isothermal condition 


Curing 


Temperature 
Tc (K) 

Cl 

C2 

f 9 

<x f (l/K) 

360 

59.0 

300 

7.36E-3 

2.45E-5 

375 

53.5 

300 

8.12E-3 

2.7 IE-5 

385 

51.8 

300 

8.38E-3 

2.79E-5 

399 

85.7 

600 

5.07E-3 

8.45E-6 

410 

110.7 

850 

3.83E-3 

4.61E-6 

425 

128.7 

1100 

3.37E-3 

3.06E-6 

435 

133.4 

1300 

3.26E-3 

2.51E-6 
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For the thermosetting resin system studied here, the magnitude of 
fractional free volume is determined by two competing factors. On one 
hand, higher curing temperatures create larger dilatational effects, 
which then generate larger free volume inside the material; on the other 
hand, increasing molecular weights as a result of polymerization 
reactions reduce the free volume due to the eliminations of number of 
molecular chain ends. The initial decrease in the fractional free 
volume fg from Tc = 360 to 385 K is, therefore, attributed to the 
increasing dilatation effects at higher curing temperatures. This is 
also consistent with the slower advancement of reactions observed in the 
chemo viscosity- time profiles in this temperature range. The decrease in 
fg from Tc = 385 to 435 K suggest that at this temperature range, the 
polymerization reaction becomes a dominant factor. 


Chapter 6 


CONCLUSIONS AND RECOMMENDATIONS 

A new analytical model for simulating chemoviscosi ty of thermo- 
setting resin has been formulated. The model is developed by modifying 
the well-established Williams-landel-Ferry (WLF) theory in polymer 
rheology for thermoplastic materials. From this investigation, the 
following conclusions and recommendation are provided: 

1. By introducing relations between the glass transition temperature 
Tg(t) and the degree of cure a(t) of the resin system under cure, 
the WLF theory is modified to account for the factor of reaction 
time. Thus the model proposed is capable of simulating viscosity 
profiles in various cure cycles. 

2. The reaction kinetics established based on the isothermal cure 
condition on the resin system is related to the case of dynamic 
curing condition. This also provides the relation for the viscosity 
in the isothermal and dynamic cases. 

3. The glass transition temperature Tg(t) in the isothermal case is 
dependent not only on the reaction time t, but also on the cure 
temperature T. Nevertheless, the relationship of glass transition 
temperature Tg ( t) and degree of cure a(t) established is also 
valid for the case of dynamic curing condition. 
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4. A single set of Cj and C 2 is found to describe adequately a material 
state which represents an average over various structures that occur 
inside during each isothermal curing condition. Thus, the ( T ) and 
C 2 ( T) may be included in the model to predict the chemoviscosity 
profiles under dynamic curing conditions. Theoretical results 
obtained from this model compare favorably well with the 
experimental data. 

5. It is recommended that the present research be extended to explore 
the application of the model to large scale processing. 
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APPENDIX A 


EVALUATION OF TOTAL HEAT OF REACTION AND DEGREE OF CURE 


The numerical analysis technique applied to determine the total 
heat of reaction and degree of cure is summarized in this appendix. 

A thermogram obtained from the DSC under the dynamic heating 
conditions is illustrated in Fig. A.l. The heat of reaction, H, can be 
expressed mathematically as 

h T - f $ dt (A.l) 

l o 


However, Hy is calculated in practice by using the following formula 


, KxAxRr 
‘t *“WxT 


(A. 2) 


where 

H Heat of reaction cal/gm 

K Calibration Constant 1/inch 

A Area under the peak of the curve inch 2 

S Chart recorder speed inch/sec 

Rr Range (Sensitivity) mcal/sec 

W Sample Weight mg 


Of these parameters, all of them are instrument parameters which can be 
set up arbitrarily except the area under the peak of the curve (Fig. 
A.l). 
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Fig. A 



X 


1 Illustration of the DSC thermogram of Hercules 3501-6 in heat 
of reaction evaluation. 
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Two numerical techniques, cubic interpolating polynomial approxi- 
mation and numerical integration were used in the calculation of the 
area under the peak of the curve obtained from the measurement of heat 
of reaction by DSC. 

For the polynomial interpolation, a set of data is obtained from 
the curve in the length scale as given below: 

xj (inch) x 0 , Xj, x 2 , ,x n 

y-( (inch) y Q , y 1# y 2 , ,y n 

In order to generalize the calculation, the method of interpolation with 
nonuniformly spaced x values was employed. Similarly, we selected the 
divided difference method instead of the Lagrangian polynomial method. 
This permits one to reuse the previous computation if it is desirable to 
add or subtract a point from the set used to construct the polynomial. 
Looking at the profile of Fig. A.l, it was presumed that the data could 
be fitted by a cubic interpolating polynomial. The 3rd-degree 
polynomial can be written as follows: 

P,(x) = a + (x-x Ja. + (x-x n )(x-x,)a, + (x-x ) (x-x. ) (x-x 9 )a, (A. 3) 

J O O 1 0 12 0X23 

where a Q , aj, a 2 and a 3 were determined by the divided difference 
method. 

The divided difference is defined as 


f[x 


0 * 




0 

0 
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f[x Q , Xy x 2 ] = 


f[x lt x 2 ] - f[x Q ,x 1 ] 


f C X q *Xj » . « « »X.. ] = 


1X2 •• • jX^ ] " f [x^ 1 « « iXj^il 


x i - x o 


also 


f (x 0 ) = ^0 


It is known that 


a o " f o 


a i * «v*i ] 


a 2 * f ^ x 0 » x i» x 2^ 


a 3 - f tVW x 3 ] 


Thus, Eq. (A. 3) can be rewritten as 


P 3 ^ “ f o + (x - x o } f t x o» x l 3 + (x ’ x 0 } ^ x_x i^ f f x 0 » x 1 * x 2 ^ 

+ (x-x Q ) (x-x 1 ) (x-x 2 ) f C x 0 » x 1 » x 2 » x 3 ^ (A. 4) 

Equation (A. 4) can be applied to every four pairs of grid points. If we 
had a set of data with (n+ 1 ) points just like x.,- and y.j, we would get n- 
1/3 cubic interpolating polynomials. 


Let f(x) = P 3 U) + Rn ( x ) 

f (4) m 

E(x) = R n (x) = (x-x Q ) (x-Xj) (x-x 2 ) (x-x 3 ) ~ ( 3 + | ) | 

It is well known that when one uses the cubic interpolating polynomial, 
the magnitude of the error term is 
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E(x) » 0(h 4 ), if lim(error) = Kh 4 

h-*o 

Since the value of h is very small, E(x) = 0(h 4 ) becomes much smaller, 
then we can approximate 


f (x) * P 3 (x) 

Thus, the value of f(x) over the entire interval [* 0 ,x n ] can be 
computed. Finally, the generalized cubic interpolating polynomial can 
be expressed as 


f i (x) = f 3i-2 + (x " x 3i-2 } fCx 3i-2’ x 31-l^ 

+ (x~x 3i _ 2 ) (x-x 3i _ 1 ) f[x 3l ._ 2 , x 3i_i» x 3i ^ 


+ (x-x 31-2 )(x-x 3i ._ 1 )(x-x 3i ) ^Cx 3i _ 2 , x 3i -_ 1 , x 3i , x 3i+J ] (A. 5) 


where i = 1, 2, 3,..., (n-l)/3 

Since every desired value of all points over the entire interval 
can be obtained by cubic interpolating polynominal mentioned above, the 
numerical integration would be applied directly upon those points. 
Before conducting the numerical integration, the coordinate transforma- 
tion is needed since the baseline of the curve obtained from the 
experiments does not remain horizontal. Using the following formula 


xj * (x.. - x Q ) cosp + (y^ - y ) sinp 


y i * “ ^ x i “ x o } sin ^ + ^ y i " y o* cos ^ (A. 6) 

-1 )-f(x ) 

i = 1, 2, 3,...,n p* tan 1 — -2- 

x n x o 


where 
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a new set of data in prime coordinates can be obtained as 




x 1 x 1 

*2 * • • • * p 



* ^2 9 * * * 


Now, the divided difference and Eq. (A. 4) in the prime coordinates can 
be obtained and a series value of f(x') can be calculated. 

In calculating the area under the peak of curve, the Simpson Rule 
[11], is employed and this is given by 


*2 5 

A = / f(x) dx * ^ {f 0 + 4 f l + V f(4) (c) x 0 < 5 < x 2 (A ‘ 7) 

x o 


where 


sx f^ (l) = O(h^) a local error 
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Equation (A. 7) can be applied on every three pairs of grid points for a 
whole set of data. Thus, for entire interval (x Q , x n ), the total area 
should be calculated by the following equation 


A ■ $ <f l-l* 4f 1 +f l + l> - - I5T 2 " 4 f<4> <51 * 0 < 5 ‘ * n (A - 8) 

where — 165 " ° (0 s 0(h 4 ) a global error 

Since the value of h selected is very small, 0(h 4 ) becomes much smaller, 
and one may approximate the area by 



ifl (f i-l + 4 f i +f i+l } 


(A. 9) 




ono on on no 
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01 ME NS ION MATER(25l*X(60)*Y(601*TABLE(60*b0)«YFINER( 100 0 ) 

01 HENS ION X FI NERtl 0001 , XPR IM C 1C001*YPRIM (1000) 

C 

C... INPUT MATERIAL PROPERTIES ... 

RE AO ( 23 . 10 ) (MATER(I>*I-1*25> 

WRITE (6*20) (MATERdl. 1-1*251 
10 FORMAT I25IAA) I 
20 F0RMAT<5X,25< AAI1 

REA0I23«A|HEIGHT.RANGE.HEAT,CHAR,H0LDT*GRTHT,ALPhA.CK,CM *CN 
WRITE(6«30I WEIGHT* RANGE »HEAT*CHAR*HOLOT*GRTHT t ALPHA. CK.CM.CN 
30 F0RMATI/5X* ‘WEIGHT -‘*F6.2*‘ I MG » • « 

• /5X, ‘RANGE »**F6.2,‘ (MCAL/SECl‘. 

• /5X,‘HcAT RATE »‘*F6.2»* IK/MINI*, 

• /5X. ‘CHART REC. -‘»F6.2*‘ CINCH/SECI*. 

• /5X, ‘HOLDING TEMP. -'*F6.2»‘ (K»‘* 

• / 5X, * GRANO TOTAL HEAT »‘,F6.2»‘ C MCAL I * ♦ 

• / 5X, ‘MAX. OEC. OF CURE -‘»F7.5,» IIOOII*, 

• /5X» ‘MOOEEL LINE SLOPE -‘ «F 7. 5, 

• /5X, *EXPONTL»M «‘*F7.5» 

• /5X, ‘EXPONTL.N «• *F7.5> 

... RE AO OATA ANO PRINT ... 

RE AO *23. •> N* M» XOF F,IBETA*(X( II»Y(I)»I-1*N) 

WRITE I 6* 100 I N.M.XOFF.IBETA, I XII I.YI Iltl-l.N) 

100 FORMAT C5X*‘ THE NUMBER OF INPUT VALUES - ‘.15. 

• /5X,‘ THE ORDER OF 01VIDED OIFFERENCEO -‘.15, 

• /5 X, • X AXIS OFF SET -‘,F9.5*‘ IINCHI*. 

• /5X* ‘ REF. INOEX ROTATION ANGLE -‘*15. 

• //14X*‘TIMEC INCHIS3X* 'HEAT GENERATION! INCH! • . 

• /( 10X*F13.6* 3X.F13.61I 

... UNIT TRANFORMATION ... 

00 50 I- 1 » N 

X( U-tX( I >*XOFF>* CHAR/60. 

YCII-Y(I IARANGE/10. 

50 CONTINUE 

BETA-ATANCYJ I BETA I— Y( II l/IXC I BETA)— XII 1 1 
WRITE(6*70I BETA 

70 F0RMATI/5X. ‘ROTATION ANGLE FOR PRIME COORD.-* »E13. 51 

... COMPUTE ANO PRINT DIVIO DIFFERENCES ... 

NAX1S-1 
202 CONTINUE 
NK-N 

CALL UTABLE(X*Y*TABLE*N.H*IFLAG*NK) 

IFIIFLAG.NE.0I GO TO 101 

WRITE 16*2001 

NM1-N-1 

00 6 I —1 ♦ NM1 

L-I 

IFf l.GT.M) L-H 

6 WRITE (6* 201 1 ( TABLE I1*J1*J-1«L) 

6 CONTINUE 

200 FORMAT </ /5X» * THE OlVlOEO DIFFERENCES TABLE ARE 

201 FORMAT ( 10X* 3E 13.5 1 

... FINER FUNCTION EVALUATION. BASE CN CUBIC INTERPOLATION POLY. 

IF CNAXIS .EO. 1 }XFINER(1I-Xtl! 

IFINAXlS.EO.ilYFINERUI-mi 
IFINAXIS.E0.2IX PR I MID -XI II 
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85 


IF CnAXIS.EO.2) YPRInC1)-Y (1) 

XI NC-.05 
JND*1 

NN-CN-1I/3 

XF-XC1) 

00 400 I-1,NN 
IND- 1*3—2 
FXO-YC IND) 

XO-XCIND) 

Xl-XC IND*1> 

X2-XC 1ND+2) 

X3-X ( I ND*3 I 

401 IF CXF*XINC.GT.X3) CO TO 400 
XF-XF^XINC 

JNO- JND*1 

IF INAXIS. EO. 2) CO TO 404 
XFINERC JND)-XF 

VF INER( JND)«FX0*TABLE!IND,1) •( XF-XO) ♦TABLE* IND*1,2 )* CXF-XO) 

* ♦ I XF— X D+TABLEClND + 2,3)*! XF-XO)* ( XF-X1 > * C XF-X 2 ) 

CO TO 401 

404 XPRIMC JNO-XF 

TP RIM JNO )-FXO-*TABLE CIND,1)* CXF-XO) ♦TABLE CIND^l.Z)* CXF-XO) 

* * ( XF-X 1) +T ABL EC IND*2«3) * C XF-XO I * C XF-X1 J * CXF-X2 I 
CO TO 401 

400 CONTINUE 

IFCNAXIS.E0.2I CO TO 503 

WRITE! 6*403) JNO.C XF INE R C I ) , 1*1 * JNO ) 

WK I TE C 6* 402 ) CYFINERCI) , 1-1, JNO) 

403 FORMAT C/5X* 'NUMBER OF FINER TIME GRID POINTS -*,15, 

* /F5X, ‘FINER TIME GRID POINTS CIN INTIAL COORO.)-', 

• / C 5X »5E1 3.5 ) ) 

402 F0RMATU/5X, 'FINER FUNCTION VALUES (IN INITIAL COORO.)-*, 

• / ! 5X ,5E13.5 ) ) 

... COORDINATE TRAFORMAT ION ... 

X1-XC1) 

Y1-YC1) 

DO 500 I - 1 , N 
Xl-XC I ) 

YI-YCI) 

XC I)-CXI-X1 )♦ CO SC BET*)* I YI-Y1)*S INC BETA) 

Y ( I ) »— ( X I— X 1) *S INC BETA) ♦CYI-Y1)*C0SC BETA) 

600 CONTINUE 
NAXIS-2 
GO TO 202 

603 WR I TE C 6 , 501 ) JNO, ! XP RIM ( I ) , I -1 , JNO ) 

WR I TE C 6, 502 ) C YPRI MC I), 1-1, JND) 

501 F0RMATC/5X, ‘NUMBER OF FINER TIME POINTS IN PRIME AXIS -*,15» 

• /5X, 'FINER TIME GRID POINTS CIN NEW COORD.)-*, 

• /CSX ,5E13.5 ) ) 

502 F0RMATC//5X, ‘FINER FUNCTION VALUES CIN NEW COORD.)-*, 

* / C 5X , 5E13. 5 ) ) 

... SIMPSON RULE INTEGRATION ... 

ITOT-JND/2 
I TOTAL- 1 TOT *2 

CALL SIHPTB C I TOTAL ,AREA1 » YPR IM.XlNC ) 

AREA-AREA1*60. 

TUT ALE- A RE A/ WEIGHT *10. 

AREA-AREA*!. /CHAR* 10. /RANGE 
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WK 1 TE ( 6 « 600 ) AREA, TOTAL! 

600 FORMAT (/ /5X»* TOTAL AREA (INCH**2. IN NEW COORO. ,•* »E 13. 5 , 
* /5X» *HEAT !NCAL)-’,E13.5) 

C 

C... SPECIAL FUNCTION EVALUATION ... 

WRITE! 6,700) 

ICOUNT-2 

702 IF(I COUNT. CE. ITOTALI GO TO 799 

CALL S IMPTB II COUNT, AREA2,VPRIN,XINC) 

AREA2-AREA2*60. 

T0TALP-AREA2/ WEIGHT* 10. 

RATIO-TOTALP/TOTALE 
FUNI-YPRIMI IC0UNT*1)/ WEIGHT 
FUN2-YPRIM! IC0UNT*1) /TOTALS/ WEIGHT 
FUN3-AL0GCFUN2) 


FUN<i*R AT 1 0 

FUN5»ALQ6(F UNA I ♦! I— ALPHA ) / AL PH A* ALOG ( l -FUN4 ) 
FUN6»T0 TAlE*CK*FUNA**CN*!1-FUM) **CN 

WRITE I 6, 701 » XPRIM I I COUNT* 1) , FUN1, FUN2 , FUN3, FUN A ,FUN5, F UN6 

IC0UNT-IC0UNT*2 

GO TO 702 


700 

0 • 

701 
799 


FORflATCZ/SX,**** SPECIAL FUNCTION EVALUATION AT SELECTED 
TINE •♦•*,/ 


/10X* ' FUNl* Y" «, 

/ 1 OX » * FUN2* Y*/H(T) *, 

/10X,*FUN3» AL0G(Y"/H1T ) ) », 

/10X * ' FUNA* H ( I )/H ( T ) 

/10X,*FUN5»AL0G(H( I l/Hl T ) )-! 1- ALPHA ) /ALPHA* ALOG ( l-Hl II/HITI) 
/10X, * FUN6*H( T)*K*(H(I)/H(T) ) **M l-H( l I / (H ( T ) ) **N» , 
//2X,*X-ININI *.3X, *FUNl ' ,7X, ' FUN2* ,8X, * FUN3* ,7X, *FUNA' , 
7X, , FUN5*,7X,*FUN6 , ,/2X,' 




FORMAT! /IX, (2F8.5I ,5! IX, El 1.5)1 

CONTINUE 

STOP 


END 


ZZZZZ22ZZZZZZZZZZZZZZZ2Z2ZZZZZZ2Z2Z2ZZZXZZZZZZZZZZZZZZZZZZZZZZZZZZ 

SUBROUTINE OTABLE! X,Y, TABLE, H,N, IFLAG.K) 

01 PENSION X (60) ,Y( 60), TABLE! 60,60) 

IF(N.LT.N) GO TO 2 
IFLAC-1 
RETURN 

2 NHl*N-l 
00 3 I-1.NN1 

3 Table! i,i)-(yu,ii-y(I))/(x( lui-xmi 

IF(N.LE.l) GO TO 6 
00 5 J-2,M 
00 5 I ■ J , NH 1 
ISU8« J ♦ 1-J 

5 TABLE! I, J)-( TABLE! I * J-l )-TABL E ( I-i, J-l) ) / (X ( I + 1 )-X ! I SUB ) ) 

6 IFLAG-0 
RETURN 
END 

CZZZZZZZZZZZZZZZZZZZZZZZZ2ZZZZZZ2ZZZ2ZZZZZZZZZZZZZZZZZZZZZZZZZZ Z2ZZZ 

SUBROUTINE SINPTB!N,A,F,0X) 

DIMENSION F ! 1 ) 




